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Abstract The rotation of erupting filaments in the solar corona is addressed through a para- 
metric simulation study of unstable, rotating flux ropes in bipolar force-free initial equi- 
librium. The Lorentz force due to the external shear- field component and the relaxation of 
tension in the twisted field are the major contributors to the rotation in this model, while 
reconnection with the ambient field is of minor importance, due to the field’s simple struc- 
ture. In the low-beta corona, the rotation is not guided by the changing orientation of the 
vertical field component’s polarity inversion line with height. The model yields strong ini- 
tial rotations which saturate in the corona and differ qualitatively from the profile of rotation 
vs. height obtained in a recent simulation of an eruption without preexisting flux rope. Both 
major mechanisms writhe the flux rope axis, converting part of the initial twist helicity, and 
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produce rotation profiles which, to a large part, are very similar within a range of shear- 
twist combinations. A difference lies in the tendency of twist-driven rotation to saturate at 
lower heights than shear-driven rotation. For parameters characteristic of the source regions 
of erupting filaments and coronal mass ejections, the shear field is found to be the domi- 
nant origin of rotations in the corona and to be required if the rotation reaches angles of 
order 90 degrees and higher; it dominates even if the twist exceeds the threshold of the 
helical kink instability. The contributions by shear and twist to the total rotation can be 
disentangled in the analysis of observations if the rotation and rise profiles are simultane- 
ously compared with model calculations. The resulting twist estimate allows one to judge 
whether the helical kink instability occurred. This is demonstrated for the erupting promi- 
nence in the “Cartwheel CME” on 9 April 2008, which has shown a rotation of ^ 1 15° up to 
a height of 1 .5 R Q above the photosphere. Out of a range of initial equilibria which include 
strongly kink-unstable (twist O = 5 tt), weakly kink-unstable (O = 3.5 tt), and kink-stable 
(0 = 2.57 t) configurations, only the evolution of the weakly kink-unstable flux rope matches 
the observations in their entirety. 

Keywords Corona, active ■ Prominences, dynamics • Coronal mass ejections, initiation 
and propagation • Magnetic fields, corona • Magnetohydrodynamics 


1. Introduction 

The geoeffectiveness of solar coronal mass ejections (CMEs) depends primarily on two 
parameters, the velocity and the magnetic orientation of the CME at the impact on the Earth’s 
magnetosphere. The higher the CME velocity and the closer its front side magnetic field 
to a southward orientation, the more intense the interaction will typically be. Therefore, 
understanding the physics that determines these CME parameters at 1 AU is one of the key 
issues in space weather research. This involves the formation and main acceleration of the 
CME in the solar corona, as well as its propagation through the interplanetary space. The 
particulars of the trigger process also play a role in some events. It appears that typically the 
corona is the place where the basic decisions are made: will the CME be fast or slow, and will 
it keep the orientation given by the source, i.e., will its magnetic axis remain oriented nearly 
parallel to the photospheric polarity inversion line (PIL), or will it rotate substantially? 

In the present paper we employ the technique of MHD simulation to carry out a first 
systematic, but in view of the problem’s complexity necessarily incomplete investigation of 
a number of processes that cause and influence changes of CME orientation in the corona. 
Such changes can be described as a rotation of the CME volume, more specifically of the 
magnetic axis of the flux rope in the CME, about the direction of ascent. This rotation should 
be distinguished from the possible rotation of the flux rope about its own axis, referred to as 
the roll effect (Martin, 2003; Panasenco et al., 201 1), which we do not address here. 

Understanding the rotation of erupting flux ropes in the corona is also relevant for the 
question which processes trigger the eruptions, as a substantial rotation may indicate the 
occurrence of the helical kink instability (KI); see, e.g., Rust and Kumar (1996), Romano, 
Contarino, and Zuccarello (2003), and Rust and LaBonte (2005). This instability is one 
of the candidate mechanisms for the initiation of CMEs (Sakurai, 1976; Fan and Gibson, 
2003; Kliem, Titov, and Torok, 2004). It commences when the twist of the rope exceeds 
a critical value, O = 2n N > <t> cr , where N is the winding number of the field lines about 
the rope’s magnetic axis. The dynamical evolution of the instability has shown very good 
quantitative agreement with a number of well observed events, which range from confined 
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filament eruptions to the fastest CME on record (Torok and Kliem, 2005 ; Williams et al. , 
2005). However, Isenberg and Forbes (2007) have pointed out an alternative mechanism for 
the rotation of line-tied flux ropes, which relies on the presence of an external toroidal field 
component, B et , due to sources external to the current in the flux rope and pointing along 
the rope, i.e., an external shear- field component. The mechanism can easily be understood 
in the simplified picture of a current loop in vacuum field. When the loop legs move out 
of their equilibrium position to a more vertical orientation, the cross product of the loop 
current with the shear-field component yields a sideways Lorentz force on the legs, which 
is antisymmetric with respect to the vertical line that passes through the apex of the loop. 
This torque forces the rising top part of the loop to rotate. The effect is also found in a full 
fluid description (Lynch et al., 2009). For a given chirality of the erupting field, it yields 
the same direction of rotation as the helical kink. Hence, a comparative study of these two 
mechanisms is required before firm conclusions about the occurrence of the KI can be drawn 
from observations of flux rope rotations, which is a further main objective of this paper. 

Since the rotations caused by the KI and by the external shear field point in the same 
direction, they are difficult to disentangle. In fact, from a more general perspective, they are 
of similar nature. Both cause a writhing of the flux rope which, by conservation of magnetic 
helicity, reduces the twist of the rope field lines about the writhing axis. Consequently, one 
can expect that observed flux rope rotations are often consistent with a range of d>-B e t 
parameter combinations which give the writhing of the flux rope by the helical kink and by 
the shear field different individual but similar combined strengths. 

Other causes of flux rope rotation include magnetic reconnection with the ambient field 
(Jacobs et al., 2009; Shiota et al., 2010; Cohen et al, 2010; Thompson, 2011; Vourlidas 
et al., 201 1), the straightening from an initial S shape (e.g., Torok, Berger, and Kliem, 2010), 
and the propagation through the overlying field. The latter comprises any asymmetric deflec- 
tion of the rising flux from radial ascent, e.g., by adjacent coronal holes (see, e.g., Panasenco 
et al., 2011), and the interaction with the heliospheric current sheet (Yurchyshyn, 2008; 
Yurchyshyn, Abramenko, and Tripathi, 2009). 

One may conjecture that the generally changing orientation of the PIL with height in the 
corona acts similarly to the heliospheric current sheet at larger heights, i.e., that the upper 
part of the rising flux continuously adjusts its orientation to align with the PIL. If this were 
the dominant effect, the rotation of erupting flux could be predicted rather straightforwardly 
from extrapolation of the photospheric field, since the overlying field is often close to the 
potential field. However, this conjecture is not valid in the lower corona where 1, and 
where the main part of the total rotation often occurs. We demonstrate this in Appendix A. 

The amount of rotation depends on the individual strengths of the five potentially con- 
tributing processes. Four of them are controlled by more than a single parameter. This is 
immediately obvious for the torque by the shear field, which must depend on the height pro- 
file B et (z), and for the reconnection, which is sensitive to the structure of the ambient field, 
i.e., whether the field is bipolar, quadrupolar, or multipolar and whether the orientation of 
the line between the resulting new footpoints of the erupting flux differs strongly from the 
original orientation. The rotation by the KI does not only depend on the initial flux rope 
twist, O — O cr , but also on the strength and height profile of the overlying field (Torok, 
Berger, and Kliem, 2010). If the overlying field decreases only slowly with height, then the 
upward expansion develops slowly and, accordingly, its contribution to the relaxation of the 
field line tension is initially weak. The relaxation is then primarily accomplished by a strong 
rotation at small heights. In the opposite case of very strong upward expansion, the rota- 
tion is distributed across a large height range, which also increases the likelihood of further 
changes by the onset of reconnection (see Lugaz et al, 2011 for an example). The effect of 
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the heliospheric current sheet can be expected to depend on the angle with the top section of 
the flux rope’s axis, on the horizontal elongation of the CME (whether its horizontal cross 
section is very elliptical or more nearly circular), and on the magnetic pressure of the CME 
relative to the pressure of the interplanetary plasma. 

Moreover, the total rotation experienced by an erupting flux rope likely depends also on 
the dynamics of its evolution. For example, a torque strongly localized at low heights, op- 
erating on a still small loop, may hurl the flux around more efficiently than a torque which 
is distributed across a large height range. As another example, in a complex (multipolar) 
coronal environment the sequence and strength of reconnection with the ambient field may 
strongly depend upon the height profiles of the rope’s angular and rise velocities caused by 
other processes, e.g., by an ideal MHD instability. The relative velocity between reconnect- 
ing flux systems controls how strongly the reconnection with the ambient field is driven. 
Hence, quantitative studies of flux rope rotation face a very high degree of complexity. 

Here we focus on two mechanisms that can cause strong rotations in the corona, the heli- 
cal kink instability and the torque exerted by an external shear-field component. By compar- 
ing a parametric study of both mechanisms in a force-free, line-tied flux rope equilibrium 
with the data of a well observed, strongly rotating erupting prominence, we demonstrate 
that their contributions can be disentangled to some degree. We also demonstrate the very 
strong influence of the ambient potential field’s height profile on the amount of rotation by 
the KI, and briefly address the influence of reconnection between the CME flux rope and the 
ambient field on the rotation. 

This investigation was stimulated by the analysis of the strong rotation in a prominence 
eruption and CME on 9 April 2008, occasionally referred to as the “Cartwheel CME”, in 
Thompson, Kliem, and Torok (2012, in the following: Paper I). Their stereoscopic recon- 
struction revealed the height-rotation profile of the erupting filament/prominence in the core 
of a CME for the first time (Thompson, Kliem, and Toeroek, 2009). This profile provides a 
strong constraint for the numerical modeling. In combination with the further observations 
of the event, it allows us to infer the major causes of the rotation and the range of source 
parameters compatible with the data. The analysis of Paper I has given the following results 
relevant for the present study. The prominence erupted from the remnants of NOAA active 
region (AR) 10989 close to the west limb and appeared as a flux rope - a single, weakly to 
moderately twisted loop - throughout the height range covered by the STEREO EUVI and 
COR1 telescopes (Howard et al, 2008), i.e., up to 4 R Q from Sun center. It rotated coun- 
terclockwise by ~ 115° up to a heliocentric height of 2.5 R Q , where the rotation leveled 
off. Two thirds of this rotation were acquired within 0.5 R Q from the photosphere. The data 
indicate a subsequent gentle backward rotation by ~ 15° in the height range up to 3.3 R Q . In 
addition, the analysis of STEREO COR2 data in Patsourakos and Vourlidas (2011) demon- 
strated that a flux rope structure is consistent also with the three-dimensional shape of the 
CME at a heliocentric distance of 13 R Q , where it had changed its orientation by a total of 
150° ± 7° from the original one, most likely by further counterclockwise rotation. At this 
stage the erupting flux was very closely aligned with the heliospheric current sheet above 
the active region. The prominence was initially accelerated mainly horizontally along the 
filament channel. This gradually turned into a radial propagation at a position ~ 98W24S 
as seen from Earth, 15° - 20° away from the original location. The prominence experienced 
most of its upward acceleration in the heliocentric height range up to ~ 2.5 R Q and reached 
a velocity of ~400 kms -1 in the COR2 field of view. At the same time, the leading edge 
of the CME accelerated to over 700 kms -1 (Landi et al, 2010). Representative images of 
the prominence from STEREO Ahead, which had the best perspective at the structure, and 
the corresponding three-dimensional reconstructions of the location of several prominence 
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Figure 1 Images and plots of the prominence eruption at 10:16 UT, as seen by the EUVI -Ahead telescope 
in the 304 A channel, and at 10:55 and 11:25 UT, as seen in white light by COR1 -Ahead. The prominence 
apex has reached heights of 0.56, 1.6, and 2.3 Rq above the photosphere at these times. The right panels 
display the reconstructed three-dimensional position of the marked prominence threads, using a reprojection 
to a viewpoint at the position of radial CME propagation, Stonyhurst longitude 98° west (relative to Earth) 
and latitude 24° south, where the counterclockwise rotation is apparent. The axes are in units of solar radii. 


threads are compiled in Figure 1 (from Paper I). The rotation (height-rotation) profile and 
the rise (time-height) profile are included below in the observation-simulation comparisons 
(Figures 6 and 8, respectively). 

As already noted above, we focus our attention here on the coronal evolution of this 
event, leaving the interaction with the heliospheric current sheet for future investigation. 
Moreover, we exclude the possible slight backward rotation by ~ 1 5° in the COR1 height 
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range from our modeling, since we are interested in the generally important effects which 
cause significant rotations in the corona. This part of the rotation is not fully certain, and, if 
real, it was likely caused by the particular structure of the large-scale coronal field above the 
active region, which nearly reversed its horizontal direction at heights > 0.3 R Q above the 
photosphere (Paper I). Thus, we will consider a saturation of the modeled rotation at angles 
near 115° and heights h ~ (1.5 -2.3) R Q above the photosphere to be in agreement with 
the observations. Furthermore, we will disregard the initial nearly horizontal motion of the 
prominence along the PIL. 

The combined effects of flux dispersal and foreshortening in the course of the source re- 
gion’s rotation to the solar limb made it impossible to obtain a well-defined estimate of the 
distance between the main flux concentrations in the bipolar region at the time of the erup- 
tion, which is a parameter of strong influence on the height profile of the ambient potential 
field. Only a relatively wide range of ~ (40- 150) Mm could be estimated by extrapolating 
the region’s evolution in the course of its disk passage through the final three days before 
the event. It will be seen that this range still sets a useful constraint on the modeling. 

In the following we model the radial propagation of the prominence in the Cartwheel 
CME in the coronal range of heights as the MHD evolution of an unstable force-free and 
line-tied flux rope (Section 2). A parametric study of the resulting rotation and rise, focus- 
ing on the rotation caused by the helical kink instability and by the external shear field, 
is compared with the observations, to constrain the parameters in the source of the event 
and to study whether the relative importance of these mechanisms can be disentangled and 
individually estimated (Section 3). The discussion in Section 4 addresses the simplifying 
assumptions made in the modeling and differences to earlier relevant work. Section 5 gives 
our conclusions. Appendix A relates the rotation of erupting flux ropes in low-beta plasma 
to the changing orientation of the PIL with height, and Appendix B considers options for 
inferring the occurrence of the helical kink in the presence of shear-field-driven rotation. 


2. Numerical Model 

We carry out a series of MHD simulations similar to the CME simulation in Torok and Kliem 
(2005). The prominence is modeled as a section of an approximately force-free toroidal cur- 
rent channel embedded in external current-free (potential) field, which represents a modifi- 
cation of the approximate force-free equilibrium by Titov and Demoulin (1999). The current 
channel creates a flux rope structure of the magnetic field which has a somewhat larger cross 
section than the channel and is enclosed by a quasi- separatrix layer in the interface to the sur- 
rounding field of arcade structure. The chirality of the flux rope is chosen to be left handed, 
so that the rotation will be counterclockwise (Green et al., 2007). The poloidal component 
of the external field, B t p , is due to a pair of subphotospheric magnetic point sources, which 
produce a pair of flux concentrations (“sunspots”) on the sides of the flux rope (the “promi- 
nence”) in the magnetogram. This field component holds the current channel in equilibrium; 
its strength at the position of the rope is exactly proportional to the current in the rope. Con- 
sequently, only its spatial profile, determined by the spacing between its sources, can be 
freely varied. The toroidal component of the external field, B et , is due to a pair of subphoto- 
spheric dipoles, positioned under the footpoints of the flux rope such that the field lines of 
B el run parallel to the magnetic axis of the rope to a very good approximation. Therefore, 
B et introduces only very minor Lorentz forces in the initial configuration, which quickly 
decrease by numerical relaxation at the beginning of each run, so that the strength of B et can 
be chosen freely within a wide range. We will also refer to the external toroidal field as the 
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Figure 2 Visualization of the modified Titov-Demoulin flux rope equilibrium used as the initial condition 
in the simulation runs of this paper; here with an average twist <f> = 3.5 jt as in Figure 4. The current channel 
is rendered as a yellow transparent volume. Blue held lines run near the magnetic axis of the flux rope (where 
the local twist is 2n ; note that their color is modified by the yellow transparent volume around them), red 
field lines are embedded in the flux surface at a distance to the axis where the local twist equals the average 
twist. Two symmetric sets of green and olive field lines show the ambient potential field. Contours of the 
magnetogram, B z (x, y, 0), are shown in the bottom plane. The torus of major radius R and minor radius a 
is submerged by a distance d, resulting in the apex height ho = R — d and the distance of each flux rope 
footpoint from the origin Df = (R- — d 2 ) 1 / 2 . A bipole, whose components are located at (±L, 0, —d), is 
the source of the external poloidal field component B ep ; see Figure 2 in Titov and Demoulin (1999) for its 
visualization. A pair of antiparallel, vertically oriented dipoles, placed under the footpoints of the flux rope at 
(0, ±Df, — 5/iq), provides the source of the external toroidal (shear) field component S e t- 


shear-field component. Here it decreases faster with height than the external toroidal field in 
the original Titov-Demoulin equilibrium. A visualization of the configuration is shown in 
Figure 2. 

We integrate the ideal MHD equations but neglect pressure, as appropriate in the active- 
region corona, and gravity, because the hydrostatic pressure profile along the field is not 
essential for the flux rope rotation, which is driven by the Lorentz force. These simpli- 
fications yield maximum freedom in the scalability of the simulation results to the data. 
Magnetic reconnection can occur due to the numerical diffusion of the field in regions of 
strong gradients. The initial density is specified as po(x) = B 0 (x) | 3//2 , where B 0 (x) is the 
initial magnetic field. This yields a slow decrease of the Alfven velocity with height, as in 
the corona. The box is a cube 64 units long on each side, significantly larger than in our 
previous simulations and in each direction at least twice as large as the biggest size of the 
structures that will be compared to the data. It is resolved by a nonuniform, fixed Cartesian 
grid with a resolution of 0.04 units in the central part of the box (a factor of two coarser than 
in Torok and Kliem, 2005). Rigid boundary conditions are implemented at the top and side 
boundaries, while very small velocities are permitted in the bottom boundary. Initially the 
torus lies in the plane {x = 0}. The MHD variables are normalized by the initial apex height 
of the flux rope axis, /z 0 , by the initial field strength B 0 , density po, and Alfven velocity 
Vao at this point, and by the corresponding quantities derived thereof, e. g. , the Alfven time 
r A = ho/ Vao- Thus, the initial apex height of the axis of the current channel and flux rope 
serves as the length unit. 

The parameters of the initial configuration are largely chosen as in Torok and Kliem 
(2005). We fix the major radius of the torus at R = 1.83, the depth of the torus center at 
d = 0.83 and the pre-normalization strength of the point sources at q = 10 14 Tm 2 in all runs. 
For a base set of the simulation series, discussed below in Figures 3 -6, 8, and 9, we further 
fix the distance of the point sources from the z axis at L = 0.83 (in units normalized such that 
ho is unity). This value lies in the middle of the estimated range for the corresponding dis- 
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tance of the flux concentrations in AR 10989, given above, when the scaling ho = 0.077 R Q 
adopted in Section 3.1 is applied. It also agrees with the settings in several previous in- 
vestigations (e.g., Torok and Kliem, 2005; Torok, Berger, and Kliem, 2010), facilitating 
comparisons. Variations of this parameter will be considered in the range L = 0.42-2.5. 
We vary the minor radius of the toroidal current channel, a = 0.32-0.62, and the strength 
of the external toroidal field, B e{ / B tp = 0- 1.06 at the flux rope apex x = (0, 0, 1), to obtain 
a range of values for the average twist of the current channel, <t> = (2.5-5.0)7t, and for 
the strength of the shear-field component. The twist is influenced by both a and B cl , with 
a having the stronger influence within the considered range of parameters. The twist values 
quoted in this paper represent the initial twist averaged over the cross section of the current 
channel in the manner described in Torok, Kliem, and Titov (2004). 

The range of the initial average twist is chosen such that unstable and stable configura- 
tions with respect to the helical kink mode are included. The first group is unstable from 
the beginning of the simulation. Nevertheless, a small upward initial velocity perturbation is 
applied in the vicinity of the flux rope apex (typically ramped up to 0.05 Vao over 5 r A and 
then switched off), to ensure that the instability displaces the apex upwards, i.e., downward 
kinking is excluded in these runs which are intended to model CMEs. 

For the geometric parameters of the system specified above, the flux rope is initially sta- 
ble with respect to the torus instability (Kliem and Torok, 2006; Torok and Kliem, 2007). 
However, the helical kink instability lifts the rope into the torus-unstable range of heights 
(h> 2 ho), from where the torus instability accelerates its top part further upwards. 1 The 
kink-stable cases require that the upward velocity perturbation is applied for a longer time, 
lifting the apex into the torus-unstable range. This allows us to study the influence of the 
shear field on the rotation in the absence of the helical kink instability, using uniform geo- 
metrical parameters of the initial flux rope (except for the minor flux rope radius a) in all 
runs. An initial velocity perturbation very close to the required minimum value is applied 
in each of these cases, to ensure nearly uniform conditions at the onset of the instability 
throughout the series. The values at the end of the ramp phase stay below 0.12 V A o for all 
runs. The flux rope velocity falls back to a much smaller value (typically ~ 0.01 V A o) imme- 
diately after the perturbation is switched off. The growing instabilities then accelerate the 
apex to peak upward velocities in the range max{« a } ~ (0.4 -0.7) V A o, far higher than the 
initial perturbation. 

On the Sun, the initial lifting of the flux can occur by a variety of effects in addition 
to the helical kink mode, as has been demonstrated by numerical simulations. These in- 
clude the shearing and twisting of the coronal field by photospheric flows (e.g., Mikic and 
Linker, 1994; Torok and Kliem, 2003), reconnection associated with flux cancellation in the 
photosphere (e.g., Aulanier et al., 2010; Amari et al, 2010), and reconnection with newly 
emerging flux (Chen and Shibata, 2000). 

The observations of the Cartwheel event indicate a gradual doubling of the prominence 
height prior to the eruption (Paper I). The initial lifting of the flux rope apex in the simula- 
tions due to the applied perturbation is much smaller for all kink-unstable runs and stays in 
the range up to this value for the kink-stable cases, except for the mn with the highest shear 
field (<J> = 2.5rr, B et /B ep = 1.06), which requires a lifting to 2.6 ho. 


1 The torus instability can be considered as a lateral kink of the current channel. However, we choose “kink” 
and “KI” to refer exclusively to the helical kink mode in this paper. 
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3. Comparison of Simulations and Observations 

3.1. Dependence of Flux Rope Rotation on Twist and Shear 

We begin with a case that involves a clear helical kink instability, as one would expect at 
first sight from the considerable rotation observed in the Cartwheel event. The initial average 
twist is chosen to be O = 5rr, a value used previously in the successful modeling of several 
filament/prominence eruptions (Williams et al, 2005; Torok and Kliem, 2005). Even with 
this considerable amount of twist (and with the sunspot semi-distance L = 0.83), we find 
that a shear field is required to reach the observed rotation. Figure 3 shows the resulting 
rotation of the flux rope, which reaches the observed value of 1 15° and is a combined effect 
of the helical kink instability and the shear field. The field lines visualize a flux bundle in the 
core of the rope which runs slightly (^ 5 %) under the rope axis in its top part. This is a likely 
location for prominence material within a flux rope. Moreover, this is the only selection that 
allows a favorable comparison with the observed flux rope shape for the weakly twisted case 
shown below in Figure 5, while the more strongly twisted cases are less sensitive to this 
vertical offset. Therefore, we adopt this selection as a uniform choice for Figures 3, 4 and 5 
which compare the flux rope rotation for different twist values. The field lines are displayed 
from perspectives identical to the STEREO images and reconstructions in Figure 1 . 

Two characteristic morphological features apparent in the COR1 data in Figure 1 are 
weakly indicated in the simulation: the initial teardrop-like appearance and the elongated 
shape at large heights (relatively narrow in the horizontal direction). The right panels show 
that the teardrop shape is a projection effect. The legs of the erupting rope approach each 
other near the edge of the occulting disk only in projection; they are displaced along the line 
of sight and actually moving away from each other. The elongated shape is largely also due 
to the strong rotation. 

The legs of the rope appear “wiggly”, which results from two effects. First, they re- 
connect with the ambient field in the vertical current sheet under the flux rope apex in the 
interval t & (32 - 65) ta, which corresponds to apex heights h ~ (5 — 21) ho', with the recon- 
nection proceeding at much lower heights inside the edge of the COR1 occulting disk. This 
leads to a bend in the reconnected flux rope: the field lines have relatively small curvature 
within the legs of the expanded original rope above the reconnection point but run along a 
more helical path in the ambient field just outside the original rope below the reconnection 
point. This bend and the more helical shape of the field lines below it relax upward, along 
with the overall upward expansion of the reconnected flux rope. Since the flux rope apex has 
reached a considerable upward velocity, n a < 0.5 Va , the bend needs a large height range 
for its propagation to the top of the rope. It is located slightly above the dotted line in the 
third snapshot pair of Figure 3 and at h > 15 ho in the final snapshot pair. The plots on the 
right hand side show that the new footpoints of the rope are displaced in counterclockwise 
direction from the original ones, thus contributing to the overall counterclockwise rotation 
of the rope. However, this contribution is only a minor one; the major part of the total ro- 
tation occurs before the flux rope legs reconnect (which can be seen by comparison with 
Figure 6 below). This reconnection is similar to the second and third reconnections de- 
scribed in Gibson and Fan (2008, their Section 4.1) and will be addressed in more detail in 
a future investigation. Second, at the given relatively high value of the twist, the dominant 
wavelength of the helical kink mode is considerably shorter than the flux rope, so that the 
characteristic helical shape develops clearly. 

Figure 4 shows the evolution in a second run where the KI develops only weakly, using a 
moderate, only slightly supercritical value of the initial average twist, <t> = 3.57T. A stronger 


Springer 


146 


B. Kliem et al. 


Figure 3 Snapshots of an 
erupting and rotating, strongly 
kink-unstable flux rope. The 
initial average twist is <t> = 5 jt 
and the shear-field component at 
the initial flux rope apex position 
is given by £ e t/#ep = 0.42. Field 
lines in the core of the rope, 
traced downward from the apex, 
are shown in the height range 
0 < z < 30, using the same two 
perspectives as for the 
observations in Figure 1 (in the 
left panels the line of sight makes 
an angle of 26° with the y axis, 
and the z axis is tilted away from 
the observer by 8° , while the 
right panels present a vertical 
view with an initial angle 
between the flux rope axis and 
the east-west direction of 26°). 
The magnetogram, B~(x,y,Q,t), 
is displayed in grayscale (seen 
from below in the left panels). 
The dotted line indicates where 
the edge of the COR1 occulting 
disk is located if the distance 
between the flux rope footpoints 
in the simulation, 2Df = 3.3 ho, 
is scaled to the value of 175 Mm 
estimated in Paper I. Using this 
scaling, the simulated heights of 
h = 1, 7.3, 21, and 30 ho (at 
t = 0, 36, 64, and 84 r^) 
translate to heights of 0.077, 

0.56, 1.6, and 2.3 Rq above the 
photosphere, reached at 10:16, 
10:55, and 11:25 UT (for rows 
2-4), respectively. 



shear field is chosen, so that the same total rotation is achieved. The overall properties - 
accelerated rise into an ejection (CME) and very strong rotation - are identical to the run 
shown in Figure 3. The morphological details, such as the teardrop shape at small heights 
and the elongated shape at large heights, match the data slightly better. The indications of 
wiggly shape at large heights remain weak. Reconnection of the flux rope legs with the 
ambient field occurs here as well, but the resulting changes in the shape of the flux rope are 
weaker, since not only the field lines in the rope are less twisted but also the ambient field 
is less helical, due to the larger B et . This morphological difference to the strongly twisted 
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Figure 4 Same as Figure 3 for a 
weakly kink-unstable case with 
initial average twist <J> = 3.5jr 
and shear field S e t/5ep = 0.67. 
The flux rope is shown at the 
simulation times t = 0, 50, 80, 
and 97 which yield the same 
heights as the snapshots in 
Figure 3, corresponding to the 
same observation times. 



flux rope is one aspect that may allow to distinguish between rotations with strong and weak 
involvement of the helical kink in observed events. The field line shapes in the present case 
conform slightly better to the inclination of the prominence threads with respect to the axis 
of the flux rope in the COR1 data in Figure 1, but this difference is not sufficiently clear 
to be decisive by itself. Moreover, it depends to a considerable degree upon which part of 
the erupting flux was outlined by prominence material in the considered event and on the 
selection of field lines in the plots. 
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Figure 5 Same as Figure 3 for a 
kink-stable case with initial 
average twist = 2.5jr and 
shear held B e t/#ep = 1.06. The 
flux rope is shown at the 
simulation times t = 0, 77, 109, 
and 128 which yield the same 
heights as the snapshots in 
Figure 3, corresponding to the 
same observation times. 



Figure 5 presents a case with subcritical flux rope twist, <t> = 2.5 it , where the kink insta- 
bility cannot develop and an even stronger shear field is needed to achieve a similar rotation. 
Here the parameters were chosen such that the rotation matches the observations as well 
as the other two runs in the height range h < 20 ho, with the total rotation of the rope’s 
magnetic axis at h = 30 ho exceeding the rotation in those runs by 20-25 degrees. The 
elongated teardrop shape at intermediate and large heights yields the best match of the three 
runs shown in Figures 3, 4 and 5. However, this is only the case because a flux bundle 
slightly under the magnetic axis of the flux rope is selected in the visualization. If instead 
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a set of field lines encircling the flux rope axis is chosen, then the high total rotation at the 
apex height h = 30 ho leads to an inverse teardrop shape (narrow at the apex, because at this 
point the view is nearly along the axis of the rotated flux rope), which is inconsistent with the 
observations. Again, since it is not known which parts of the erupting flux (rope) were filled 
with prominence material in the event to be modeled, these morphological comparisons, by 
themselves, do not allow to rule out the kink-stable run shown in Figure 5. 

The similar total rotations in the three simulations confirm that both twist and shear 
belong to the key parameters which determine the amount of rotation in erupting flux ropes. 
To analyze this further, we consider a set of characteristic cases from our series of simulation 
runs with varying strength of the two effects. For each of the twist values O = 5.0, 3.5, and 
2.5jt , we vary the shear field B cl from the respective best fitting value used in Figures 3-5. 
All runs use the same sunspot semi-distance L = 0.83 and, hence, the same external poloidal 
field B ep . The variation of L will be considered in Section 3.2. 

The rotation of the flux rope in the simulations is measured in two ways. At low heights 
it is taken from the changing orientation of the magnetic axis at the apex of the flux rope. 
As the flux rope rises, the apex orientation oscillates increasingly, due to the upward propa- 
gation of Alfvenic perturbations which result from the dynamic onset of reconnection in the 
vertical current sheet under the rope (the relaxation of the bend in the reconnected field lines 
mentioned above). The right panels at the two final heights in Figures 4 and 5 indicate the 
resulting oscillations of the field orientation at the apex with respect to the bulk orientation 
of the flux rope’s upper part. Therefore, at larger heights we simply use the direction of the 
horizontal line connecting the flux rope legs at the height where they are most distant from 
each other. This measurement filters away most of the oscillating variations, which are also 
not captured by the observed rotation data derived in Paper I and replotted in Figure 6. The 
difference between the two measurements remains less than five percent in a height range 
Ah ~ (3 - 6) ho around h ~ 10 ho, except for the most strongly rotating and oscillating case 
in the series (<f> = 57 t, B et /B ep = 0.63) where it reaches ~ 10 percent. Linear interpolation 
between the two measurements for each simulation run is applied in the appropriate range 
of small differences to match them smoothly. 

(The method to estimate the rotation angle at large heights fails for one of the runs in 
Figure 6(0 = 2.5jt, B et = 0), where reconnection of the flux rope legs with the ambient 
field leads to jumps that are larger than the oscillations of the magnetic axis at the apex. 
For this run, whose rotation profile differs strongly from the observed one, we include the 
rotation angle only at low heights, to show the trend.) 

In order to compare the simulated rotation profiles with the observations, a scaling of the 
length unit in the simulations to distances on the Sun is required. For this purpose, we set the 
distance between the footpoints of the flux rope in the simulation, 2D{ = 3.3 ho, equal to the 
estimated length of the flux which holds the prominence, 175 Mm (Paper I). This is indepen- 
dent of the actual prominence shape. The apex height of the toroidal Titov-Demoulin flux 
rope, our length unit, tends to be somewhat high in comparison to solar prominences, which 
are often quite flat. Here we obtain ho = 0.077 R Q , relatively close to the estimated initial 
prominence height of ~ (0.05-0.06) R Q (Paper I). If we would instead choose to compare 
the simulations to the temporal profile of the prominence rotation, then each change of the 
twist, which implies a change of the KI growth rate, would require a rescaling of the time 
unit in the simulations, The comparison of the simulated rotation profiles with the ob- 
served profile is displayed in Figure 6. As discussed in Sections 1 and 4, we disregard the 
slight backward rotation at h > 1.5 R Q above the photosphere in the comparison and as- 
sume that the tendency of the rotation to level off at this height would have continued in the 
absence of the specific complex structure of the large-scale coronal field above AR 10989 
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Figure 6 Comparison of flux rope rotation as a function of normalized apex height above the photosphere 
with the observations obtained in Paper I. Crosses and diamonds are EUVI data, with the final cross rep- 
resenting a lower limit for the height and the diamonds representing interpolated heights. Plus symbols are 
COR1 data. The distance between the footpoints of the flux rope in the simulation is scaled to the value of 
175 Mm estimated in Paper I, resulting in h q = 0.077 Rq . The initial average twist, O, and the strength of the 
shear- field component (external toroidal field), B e t, given by its ratio to the external poloidal field component 
B e p at the initial flux rope apex, are varied, while the geometrical parameters of the initial flux rope (except 
the minor radius a) and the spatial structure of the external field components B e t and B e p are uniformly cho- 
sen throughout the series of runs (see Section 2 for details). The optimum values for the shear-field strength, 
which yield the best match with the observed rotation profile up to h » 20 /zq, found through parametric 
search, are B e t, opt/^ep = 0.42, 0.67, and 1.06 for <1> = 5.Cbr, 3.5 jt, and 2.5 j r, respectively. Changes of B e t 
by a factor 3/2 and the case B e t = 0 are included. 


and in the absence of the heliospheric current sheet, which are not included in our model. 

Several conclusions can be drawn from this set of simulations. 

i) Similar height-rotation profiles (not only a similar total rotation) are obtained in a range 
of O-tSet combinations. The profiles for (d>, B et /B e p ) = (57r,0.42), (3.5^, 0.67), and 
(2.571, 1.06) all match the observed profile very well up to a height h ~ 20 ho ~ 1 .5 R Q 
above the photosphere, where a total rotation of ^ 1 15° is observed. These runs include 
a strongly and a weakly kink-unstable and a kink-stable case. Hence, even such a strong 
rotation does not by itself imply the occurrence of the helical kink instability. Further 
arguments, such as those given below, are required to draw conclusions about the oc- 
currence of the instability in the modeled event. 

ii) To reach the observed total rotation of & 1 15° with the initial configuration and parame- 
ter settings chosen in this series, in particular with the chosen value of the sunspot semi- 
distance L, the shear must contribute. The strongly twisted configuration (O = 5jt) 
yields only little more than one third of the observed rotation in the absence of shear 
( B et = 0). Therefore, the shear contributes the main part of the total rotation even in this 
strongly kink-unstable case. Note that this conclusion changes if the sunspot distance 
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is set to larger (however, unrealistic) values, so that the overlying field decreases less 
steeply with height (see Section 3.2). 

iii) The twist also contributes in all runs. The tension of the twisted field relaxes in any case 
when the flux rope is driven upward out of its initial equilibrium, be it by the helical 
kink instability, by the torus instability, or by any other process (e.g., by so-called tether- 
cutting reconnection). This relaxation contributes to the writhing of the flux rope axis 
regardless of whether or not the helical kink instability is triggered. As a consequence, 
we do not observe a jump in the achieved rotation as the twist of the initial equilibrium 
is varied between kink-stable and kink-unstable values. This is most obvious from the 
runs with B et = 0. 

iv) The higher the relative contribution of the twist, the lower the height range where most 
of the rotation is reached. This reflects the fact that the KI tends to reach saturation 
quickly, often already when the flux rope has risen to a height comparable to the foot- 
point distance (e.g., Torok, Kliem, and Titov, 2004). This property corresponds well to 
the tendency of the rotation to level off at the relatively low height of ^ 1 .5 R Q (~ 20 ho) 
above the photosphere. The rotation by the shear field acts in a larger height range. The 
different behavior can be made plausible from the fact that the Lorentz force due to the 
shear field depends on the current through the rope and on the angle between the flux 
rope legs and the shear field. While the current decreases as the rope ascends (similar to 
the twist), the angle rises until the legs approach a vertical position, which corresponds 
to bigger apex heights than the saturation height of the helical kink mode. Hence, the 
Lorentz force due to the shear field acts strongly in a larger height range than the tension 
force associated with the twist. 

As a consequence, the Titov-Demoulin flux rope with subcritical twist for KI onset 
does not allow to match the entire observed rotation profile of the 9 April 2008 event. 
We have performed considerable numerical experimenting in this range of twists [O = 
(2.5-3)7t], including modifications of the height profiles B et (z) and B ep (z) and of the 
flux rope shape (by varying its major radius R but not the apex height ho) from the 
uniform settings for the runs in Figure 6. Either the rotation in the height range h < 20 ho 
was found to be too small, or the total rotation at A = 30 ho was too large. Although the 
shape of the prominence in the plane of the sky can still be met by the special selection 
of the field lines in Figure 5, the saturation of the rotation at A 1.5 2? 0 , revealed by 
the stereoscopic reconstruction, cannot be reproduced. This suggests that at least a weak 
helical kink instability must have been triggered in this event. 

v) The range of twist-shear combinations that reproduce the observed rotation profile is 
bounded not only from below, as outlined in ii) and iv), but also from above. Average 
twists significantly exceeding 5 it are not only unlikely to occur in the corona but also 
lead to increasingly strong helical deformations of the flux rope, which are favorable 
for the onset of magnetic reconnection with the overlying field or between the flux rope 
legs. Such reconnection can strongly distort the rotation profile and can even stop the 
rise of the flux rope (Torok and Kliem, 2005 ; Shiota et al. ,2010). Reconnection with the 
overlying field does indeed lead to a confined (failed) eruption in the present simulation 
series when the initial twist is raised to 6rr . Reconnection between the legs of the rope 
occurs if O > Ire , also leading to confined eruptions. (A detailed description of such 
reconnection can be found in Kliem et al., 2010.) 

Increasing the shear field tends to stabilize the flux rope because any displace- 
ment then requires an increasing amount of energy to push the ambient field aside. 
The low-twist case (O = 2.5tt) with the strongest shear field included in Figure 6 re- 
quires a considerable initial perturbation to reach the torus-unstable range of heights 
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(. h > 2.6 ho = 0.2 R q for these parameters); it is completely stable to small perturba- 
tions. Similarly, while the 3.57T run with B a = 0 is clearly kink-unstable, the corre- 
sponding sheared case (B et /B ep = 0.67) exceeds the instability threshold only slightly. 
The initial lifting of the flux rope required in the low-twist case strongly exceeds the 
observed rise of the prominence to ~ 0.06 R Q prior to the onset of the eruption. This 
represents a further strong indication against this configuration. 

The upper limit for the shear field is not a universal number but depends on other 
parameters of the system, which include the thickness of the flux rope, the strength of 
the line tying, and the height profile of the external poloidal field, 5 ep (z). A systematic 
study of these dependencies would be beyond the scope of the present investigation. 
However, we have considered a modification of the height profile B ep (z), which is the 
key parameter for the onset of the torus instability in the absence of shear and significant 
line tying (Kliem and Torok, 2006). In an attempt to ease the occurrence of the insta- 
bility in the low-twist case (<t> = 2.5rr, B et /B ep = 1.06), the sunspot semi-distance was 
reduced to the minimum value of the possible range estimated from the observations, 
L = 0.4, leaving the other parameters of the equilibrium unchanged. No reduction of the 
minimum height for instability was found, which must be due to the strong stabilizing 
effect by the chosen shear field. 

vi) Reconnection of the flux rope legs with the ambient field contributes only a minor part 
of the total rotation in our simulation series. It appears to remain weaker than the twist- 
driven rotation, or at most comparable, i.e., considerably weaker than the shear-driven 
rotation. This can be seen most clearly in the 57 t run with B ct = 0. Here the reconnection 
of the flux rope legs with the ambient field proceeds while the rope apex rises from 
2 ho to « 16 ho, with the flux in the core of the rope being involved in the range of 
apex heights h ~ (4- 16) ho- However, the major part of the total rotation of ~ 40° is 
already reached at low apex heights, h < 5 ho, i.e., due to the helical kink mode. The 
apex height range during the reconnection of the flux rope legs in the shear-free 3.5 tt 
run is similar to the 5jt run. The rotation profile of this run in Figure 6 shows about equal 
amounts of rotation in the height ranges h <5 ho and /z ~ ( 5-16) ho, indicating that 
the reconnection-driven rotation could here be comparable to the twist-driven rotation. 
Again, both remain considerably smaller than the rotation due to the shear in the 3.5 tt 
run that best fits the observations. 

These conclusions are also supported by the fact that the angular distance between 
the initial and new footpoints of the flux rope’s magnetic axis, measured from x = 0, 
remains far smaller than the total rotation of the rope (see the right panels in Figures 3 - 

5). 

3.2. Influence of the External Poloidal Field 

The height profile of the poloidal field which is due to sources external to the flux rope, B ep , 
is a further factor of potentially strong influence on the rotation. Erupting flux ropes rotate 
more strongly at low heights if the external field initially overlying the flux rope decreases 
more gradually with increasing height (Torok, Berger, and Kliem, 2010). The relaxation of 
the magnetic tension in the erupting flux rope by rotation is then more pronounced because 
the relaxation by upward expansion is hindered, at least initially. The relevant length scale, 
l z = — [d(log fi e p)/dz] -1 , increases with increasing distance between the sources of B ep , i.e., 
between the main flux concentrations to the sides of the PIL. This can easily be seen for the 
Titov-Demoulin equilibrium, where this scale height is l z = (z + d)\ I + L 1 / (z + d) 2 ]/ 3. 

Figure 7 shows that this effect remains weak as long as the distance between the sources 
of B ep , 2 L, is smaller than the distance between the footpoints of the erupting flux rope, 
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Figure 7 Dependence of flux rope rotation vs. height upon the distance L of the main flux concentrations 
in the source region from the PIL. Kink-unstable flux ropes (<t> = 5jr and 3.5 jt) are considered for vanishing 
external shear-field component, B e t = 0. 


2Df, but that it becomes very strong when the reverse relation holds. Here the sunspot semi- 
distance L is varied for the 57 t and 3.5rr runs with no external shear field, B el = 0, to be 
0.5, 1, 2, and 3 times the value estimated from the observations and used in Section 3.1 
(Figures 3-6). The two distances are nearly equal if L is set to twice the estimated value. 
This is larger than the maximum of the range for L compatible with the observations (see 
the Introduction). Hence, the conclusions drawn from the series of simulations shown in 
Figure 6 are not sensitive to the actual value of the parameter L as long as it remains within 
this range. In particular, an external shear-field component of strength close to the optimum 
values given in this figure is then required to reach the observed rotation. 

Rotations even exceeding those produced mainly by the shear field in Figure 6 are 
achieved in the absence of a shear field for both twists if L exceeds Df by a factor >1.5. 
A similar situation was realized in simulations of erupting flux ropes in Fan and Gibson 
(2003) and Gibson and Fan (2008), which showed strong rotations of 115-120 degrees 
with B et = 0. However, such large distances of the main polarities, relative to the length 
of the PIL and a filament channel between them, do not typically occur in fully developed 
active regions. Hence, the effect of a shear field (Isenberg and Forbes, 2007) will typically 
be involved (and significant) if erupting flux rotates by large angles of order 90° and more. 

3.3. Rise Profile 

The results of Sections 3. 1 and 3.2 lead to the question whether the initial twist and the shear 
field in the source volume of the eruption can be further constrained individually, although 
their combined effect on the rotation is similar. The rotation profile obviously is a powerful 
new diagnostic of the evolution of flux ropes in CMEs, however, for the considered event it 
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does not allow to discriminate between the strongly and weakly kink-unstable cases shown 
in Figures 3 and 4, respectively. Therefore, we now consider the rise (time-height) profile of 
the erupting flux. This function reflects the growth rate of the instability driving the eruption. 
The growth rate varies strongly with the twist if this parameter exceeds the threshold of the 
helical kink mode (see, e.g., Figure 5 in Torok, Kliem, and Titov, 2004). When the variation 
of the twist is combined with a variation of the shear-field strength in the opposite direction 
(one increasing, the other decreasing), such that the rotation profile stays nearly unchanged, 
then the rise profile will change even stronger: decreasing (increasing) shear-field strength 
leads to higher (lower) KI growth rate. Thus, the combined comparison can constrain these 
parameters individually. 

In order to compare the simulated rise profiles with the observed one, the time unit in 
the simulations, r A , must also be scaled to a dimensional value. Since r A = ho/ V A o and ho 
is already scaled, this is equivalent to adopting a value for the initial Alfven velocity F A0 
in the body of the prominence. So far, this parameter can hardly be derived from observa- 
tions, since both the field and density structure of prominences are generally only poorly 
known. Therefore, here we work backwards by first finding the best match between the sim- 
ulated and observed rise profiles and then checking whether the implied Alfven velocity 
falls within an acceptable range. Lower bounds on the Alfven velocity in filaments have 
been obtained through the application of seismological techniques to six cases of oscillating 
filament threads (Terradas et al., 2008). Five of these lie in the range ~ (300-600) kms -1 
if the length of the field lines that pass through the threads is assumed to be ~ 175 Mm, the 
length of the erupting structure estimated in Paper I. An upper bound of order 1000 km s -1 
is widely accepted for old, dispersed active regions like the one considered here. 

The rise profiles of the simulation runs shown in Figures 3, 4 and 5 are scaled and 
matched to the observed profile in Figure 8. In selecting the scaling parameters for the 
best match, we adopt a start time of the eruption a couple of minutes before 08:51 UT, as 
estimated in Paper I. The conclusions drawn from the comparison do not depend upon the 
particular start time if chosen in this range. The value 08:48 UT used in Figure 8 yields the 
best match of the 3.57T and 2.5 tt runs with the observations and lies very close to (30 sec 
before) the last EUVI image prior to the occurrence of motions in the prominence along the 
path of the CME. Also, we give relatively low priority to the EUVI height data after 10 UT, 
since these may be smaller than the true heights, as discussed in Paper I. 

The scaled rise profile of the simulation with O = 3.5 tt is found to fit the data quite well 
if the Alfven velocity is chosen in the range V A0 = (540-560) kms -1 and the start time of 
the simulation is placed in the range 08:45-08:50 UT (with the earlier time corresponding 
to the lower V A o). These values appear very plausible. 

We did not succeed to find a satisfactory fit by the higher twisted case. The corresponding 
curve in Figure 8 demonstrates this, using the same start time as for the 3.5rr run and V A n = 
420 kms -1 . Increasing (decreasing) V A0 leads to a steeper (flatter) fit curve, i.e., to abetter fit 
at the larger (smaller) heights (if the start time is adjusted simultaneously), but it is obvious 
that the curve can never fit the combined EUVI and COR1 time-height data. Here the phase 
of accelerated rise ends too early because the instability grows and saturates too quickly. 
The rise profile of this simulation can be stretched on the time axis and formally be fit to the 
data if in addition to an unrealistically low Alfven velocity of 300 kms -1 (lower than the 
terminal speed of the CME core) an unrealistically large extension of the prominence flux 
of 360 Mm (twice as large as the estimate in Paper I) are assumed. Both are unacceptable. 
This comparison with the data thus argues clearly against the occurrence of high twist and a 
strong helical kink instability in the considered event, in spite of the high total rotation. 

Assuming the same start time as for the other two runs, the kink-stable low-twist case 
(<t> = 2.5 tt) allows an acceptable approximation of the observed rise profile, which yields 
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Figure 8 Comparison of the 
observed and simulated rise 
profiles of the flux rope apex, 
using the same scaling of lengths 
in the simulations as in Figure 6 
and a start time of the eruption at 
08:48:00 UT. EUVI and COR1 
data from Paper I are plotted 
using the same symbols as in 
Figure 6. The 5 it, 3.5n, and 
2.5 tt runs of Figures 3 - 5 are 
scaled to these data assuming 
Alfven velocities 
V A0 = 420 kms -1 , 550 kms -1 , 
and 560 kms -1 , respectively. 



Time (UT on 9-Apr-2008) 


a plausible value of 560 kms -1 for the Alfven velocity. The match is slightly worse in 
comparison to the 3.5 jt run because the curve does not reach the height of the first COR1 
data point. Reducing Vao> and adjusting the start time, allows for a nearly perfect match of 
the COR1 data, similar to the 3.5 jt run, but this moves the simulation curve, which already 
runs above all EUVI data points, further away from the measurements in this height range, 
so that the overall match is degraded. 

The origin of the difference lies in the tendency of the torus instability to spread the main 
upward acceleration of the flux across a larger height range than the helical kink instability, 
which can be clearly seen in Figure 8. The height range for the torus instability is small 
only if the field in the source volume of the eruption decreases very rapidly with distance 
from the flux rope position (see Figure 1 in Kliem and Torok, 2006), i.e., in very compact 
active regions of high field strength, especially in quadrupolar ones. Since AR 10989 was 
already rather diffuse by the time of the eruption, there is no justification to make the initial 
configuration in the simulations more compact for a better fit of the rise profile by the kink- 
stable configuration. 

3.4. Implications for the 9 April 2008 Eruption 

Based on the good quantitative agreement of the simulated rotation and rise profiles with 
the observations, Sections 3.1, 3.2 and 3.3 yield the following picture. The rotation profile 
in the height range h < 20 ho & 1 .5 R Q above the photosphere is well matched by a strongly 
kink-unstable case (O = 57r), a weakly kink-unstable case (<b = 3.5n), and a kink-stable 
case (<f> = 2.57 t) if a shear field of appropriate strength is included in each of them. At 
greater heights, h ~ (20-30) ho ~ (1.5 -2.3) R 0 , the comparison yields a clear indication 
against the kink-stable case, which enters this range with an accelerated rotation, while the 
observed rotation levels off. The kink-stable case also requires a considerably stronger initial 
perturbation, lifting the flux rope apex into the torus-unstable range of heights, i.e., to h > 
2.6 ho = 0.2 R q , a value not supported by the observations. In comparison, the accelerated 
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rise of the kink-unstable cases in our simulation series starts essentially from ho = 0.077 R Q , 
relatively close to the observed onset height of (0.05-0.06) R 0 . The shear field required by 
the kink-stable case is comparable to the external poloidal field, B et / B ep = 1 .06. In a bipolar 
region, this corresponds to a similar distance between the main polarities along and across 
the PIL, which is not supported by AR 10989 as long as its magnetic structure could be 
discerned in the approach to the limb (see Figure 4 in Paper I). The rise profile rules out the 
strongly kink-unstable case and yields a further indication against the kink-stable case, albeit 
only a weak one. Both the observed shape of the flux rope as a whole and the observed angles 
between individual threads and the rope axis can be approximately reproduced by all three 
model systems, but the overall match is best for the weakly kink- unstable case (Figures 3, 
4 and 5). This is substantiated by Figure 9, where we plot the sets of field lines for this and 
for the kink-stable case which were found to match the observations closest, out of many 
different sets that were considered. 

The shape of the erupting flux rope’s magnetic axis in the considered event is not suf- 
ficiently well defined by the observations to allow a clear discrimination between the three 
considered cases based on this property alone. Note that for other events it has proven to be 
decisive. For example, the shape of the two erupting filaments modeled in Torok and Kliem 
(2005) could be matched only if an initial average twist of 57 t was assumed, not with a twist 
of 4 ji. 

Overall, we conclude that both strongly kink-unstable and kink-stable configurations can 
be excluded with a high degree of certainty, leaving a weakly kink-unstable initial configura- 
tion as the most likely source of the Cartwheel event. This configuration allows to reproduce 
the event with observationally supported values for several key parameters (flux rope length, 
distance of the main flux concentrations, initial orientation) and with plausible assumptions 
for the magnetic structure (flux rope in a simple bipolar active region) and for the remaining 
free parameters (twist and shear-field strength). 

Regardless of how definite the rejection of the other two cases is considered to be, the 
rotation of the erupting flux was primarily caused by a shear field (Isenberg and Forbes, 
2007). Weaker contributions came from the relaxation of twist (most likely by a weak helical 
kink instability) and from reconnection with the ambient field. 


4. Discussion 

The major simplifying assumptions adopted for the modeling in this paper include 

i) the neglect of the initial mainly axial propagation of the prominence, 

ii) the neglect of any asymmetry and complexity introduced by the large-scale overlying 
field, and 

iii) the assumption of a well-defined, coherent flux rope (i.e., the Titov-Demoulin model). 
We discuss these here to assess their potential influence on the results. We also compare 
our results with other recent relevant work. 

While the initial propagation of the prominence introduced an asymmetry and, therefore, 
definitely had the potential to produce some rotation, we expect that it could not contribute 
strongly because the propagation was approximately along the flux holding the prominence. 
This does not principally change the magnetic configuration and the Lorentz forces which 
dominate the acceleration of plasma in the low-beta corona. 

The effects belonging to category ii) are likely to be relevant primarily at considerable 
heights. AR 10989 was a relatively isolated region of simple, bipolar structure, and this 
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Figure 9 Comparison of simulated and observed flux rope shapes for the kink-stable run (<J> = 2.5 jt, 
Bet/Sep = 1-06; left panels) and the weakly kink-unstable run (<t> = 3.5 jt, S e t/5 e p = 0.67; right panels) in 
our parametric search which best match the observed rotation and rise profiles in their entirety. The STEREO 
images from Figure 1 are supplemented by an additional image at 10:26 UT from Paper I. For both runs, 
some experimenting with the field line selection was performed until also the observed shape was matched 
best. This yielded a flux bundle running slightly under the apex point of the rope’s magnetic axis for the 
kink-stable run, as in Figure 5, and a flux bundle enclosing the axis for the kink-unstable run. 


holds also for its dispersed phase as long as it could be followed in the approach to the limb. 
The potential-field source-surface extrapolation of the photospheric field in Paper I shows 
that the large-scale coronal field associated with the polar fields and the heliospheric current 


<£) Springer 





158 


B. Kliem et al. 


sheet began to dominate already at heights h > 0.3 R Q above the photosphere, where the 
horizontal field direction nearly reversed. The force by the field component along the line 
between the flux rope legs pointed in the direction of a clockwise rotation above this height, 
opposite to the force low in the corona. However, the shear field above ~ 0.3 i? Q was weaker 
than the shear field in the core of the active region by more than an order of magnitude, so 
that it could efficiently counteract the continuing, oppositely directed force by the shear field 
at low heights, and the angular momentum of the already rotating flux rope, only by acting 
across a considerably larger height range. This is consistent with the fact that the possible 
weak reverse rotation occurred only at h > 1 .5 R Q above the photosphere. Thus, the rotation 
caused by the shear field and twist inside the bipolar active region (at h < 0.3 R Q ) must 
have been dominant factors for the rotation in the height range up to ~ 1.5 R Q modeled 
here. We cannot exclude that the saturation of the rotation would have occurred at a greater 
height if the horizontal field had not changed its direction above the active region, however, 
this weakens only one of the three main arguments against the kink-stable configuration 
summarized in Section 3.4. The saturation of the rotation profile, at a very similar height, 
was also seen in another erupting quiescent filament (Bemporad, Mierla, and Tripathi, 2011; 
see their Figure 5). 

The effect of the heliospheric current sheet is expected to become important only at even 
larger heights. Otherwise, the rotation would not have shown the saturation near h~ 1.5 i? Q 
and the possible subsequent slight reverse rotation; rather the continuation of the rotation 
to the value of ~ 150° found at 1 3 R Q would have proceeded already in the COR1 height 
range. 

The assumption that erupting flux in CMEs takes the structure of a flux rope is strongly 
supported by all available observations. Quantitative differences to our modeling must occur 
when initial flux ropes of different structure are used. These are not likely to be substantial 
if only details of the structure differ. The helical kink mode is known to not overly depend 
on the details of the current channel’s radial structure. This can be seen, for example, from 
the similar instability thresholds found in Mikic, Schnack, and van Hoven (1990), Baty and 
Heyvaerts (1996), Torok, Kliem, and Titov (2004), and Fan and Gibson (2003) although flux 
ropes with and without a net current and with straight and arched geometries were investi- 
gated. Flux ropes with hollow current channels have recently been found to be representative 
of filament channels which have undergone substantial amounts of flux cancellation (e.g., 
Su et al., 2011). It is conceivable that their less compact current distribution leads to smaller 
rotations than the Titov-Demoulin equilibrium with the same twist. This will be a subject 
of future study. On the other hand, we believe that a strongly kink-unstable configuration of 
this type would likely still not match the observed rise profile. The structure and strength 
of the external poloidal and toroidal field components do not depend upon the details of the 
flux rope structure, so that two arguments against the kink-stable configuration, which are 
based on the required initial lifting and on the ratio of B el and B ep , would likely still apply. 

An overlying current sheet (Birn, Forbes, and Schindler, 2003) may be of stronger influ- 
ence, but we have argued above that this was not the case for the considered event at the low 
coronal heights modeled in this paper. 

The situation likely changes if the flux rope is far less coherent than the Titov-Demoulin 
configuration (Green, Kliem, and Wallace, 201 1), especially if it is split (Bobra, van Balle- 
gooijen, and DeLuca, 2008). The investigation how such complex cases might change our 
conclusions must be left for future work. 

The comparison of the flux rope rotations found in this paper with the rotation in the 
simulation of a breakout CME by Lynch et al. (2009) suggests a strong dependence upon the 
existence of a flux rope at the onset of the eruption. In that simulation, the inflating flux of a 
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continuously sheared arcade did not show any significant rotation up to a heliocentric height 
of ^2 R q . Flare reconnection commenced at this point, which progressively transformed 
the inner part of the arcade into a growing flux rope. The flux rope immediately began to 
rotate. This process was monitored until the core of the rope reached a heliocentric height 
of ^3.5 Rq. Throughout this range, the rope showed a linear increase of its rotation angle 
with height, and the twist in the rope stayed below the threshold of the helical kink mode. 
The addition of poloidal flux by flare reconnection was largely complete in the middle of the 
height interval. The rotation profile in this model differs principally from the data presented 
here, even if only the height range > 2 R Q is considered, where a flux rope did exist. This 
suggests that the presence of a flux rope at the onset of the eruption was a key feature of the 
Cartwheel event. 

An interesting result of our parametric study is that the erupting flux rope did always 
show some amount of rotation, even in the shear-free, kink-stable case included in Figure 6. 
We expect this to be generally valid if coherent force-free flux ropes are considered as the 
initial condition, because such ropes always possess twist. An untwisted flux tube, known as 
a Theta pinch, requires a radial pressure gradient to attain equilibrium. This is not available 
if the plasma beta is very small, as expected for the lower coronal part of active regions. 
Whether the observations support the occurrence of rotation in essentially all events does 
not yet seem to be clear. For example, Muglach, Wang, and Kliem (2009) report that only 
about 10 cases of unambiguous rotation in erupting filaments not very far from Sun center 
could be identified in the EUV observations by the EIT instrument (Delaboudiniere et al . , 
1995) for the whole Solar Cycle 23. Flowever, many cases of only moderate rotation may 
remain undetected in such data, due to the projection on the plane of the sky. Thernisien, 
Floward, and Vourlidas (2006) find that only 13 out of 34 events fitted by a geometric flux 
rope model showed rotations (between 5° and 90°) at distances up to 30 Rq, but the ro- 
tation angles obtained through this method appear to possess a considerable uncertainty 
(Thernisien, Vourlidas, and Howard, 2009). Yurchyshyn, Abramenko, and Tripathi (2009) 
report 101 partial and full halo CMEs which show a very broad distribution of the differ- 
ence between the estimated initial and final orientations in the same distance range; these 
angles do not show a clustering at zero degrees. However, the latter two investigations yield 
the net effect of rotation in the corona and in the inner solar wind where the heliospheric 
current sheet likely dominates. If the fraction of non-rotating events is relatively small, then 
a plausible explanation is that other processes counteract the rotation by twist relaxation and 
the shear field in these cases, for example reconnection with the ambient field. If the frac- 
tion is large, then such nearly exact cancellation of rotations is unlikely to be the primary 
explanation. The implication would then be that the current distribution in the erupting field 
is often less compact or less coherent than in the Titov-Demoulin flux rope, including the 
possibility that a flux rope does not yet exist at the onset of the eruption. 


5. Conclusions 

Our parametric study of force-free flux ropes which erupt from simple bipolar source regions 
with no overlying current sheet and rotate about the direction of ascent yields the following 
conclusions. 

i) Both the force by an external shear-field component B et (Isenberg and Forbes, 2007) 
and the relaxation of twist <f> (e.g., Torok, Berger, and Kliem, 2010), are potentially 
very significant contributors to the rotation. 
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ii) For parameters typical of CME source regions, in particular if the sources of the ex- 
ternal stabilizing field (usually the main flux concentrations next to the PIL) have a 
smaller distance than the footpoints of the erupting flux, the shear field yields the dom- 
inant contribution to the rotation for a wide range of shear-field strengths. The relax- 
ation of twist remains the weaker contributor under these conditions, even if the twist 
is sufficiently high to trigger the helical kink instability. However, since twist always 
exists in force-free flux ropes, it always causes at least some rotation. Strong rotations 
(>90°) can be produced by the twist alone, but only for considerably larger distances 
between the sources of the external stabilizing field than typically observed. 

iii) The rotation in low-beta plasma is not guided by the changing orientation of the PIL 
with height. For the geometrical conditions typical of CME source regions, it is oppo- 
site in direction (see Appendix A). 

iv) For a given chirality of the configuration, the external shear field and the twist cause 
flux rope rotation in the same direction, which is clockwise for right-handed field and 
counterclockwise for left-handed field if seen from above. 

v) The two processes are related to each other when considered in terms of magnetic 
helicity. Both convert initial twist helicity of the flux rope into writhe helicity. The 
same total rotation, and rotation profiles which are very similar in a substantial part of 
the total height range of rotation, result in a range of B et -< f> combinations. 

vi) The rotation due to twist relaxation tends to act mainly low in the corona, in a height 
range up to only a few times the distance between the footpoints of the erupting flux. 
The rotation by the shear field tends to be distributed across a larger height range. 

vii) The mere fact that erupting flux rotates does not by itself imply that the helical kink 
instability occurred. Inferring this instability requires deriving a supercritical twist or 
finding features incompatible with rotation driven (nearly) exclusively by the shear 
field, like more than one helical turn of the erupting structure, approaching legs, or a 
rotation exceeding « 130° at coronal heights (in the absence of other contributions). 
A somewhat smaller rotation angle may be sufficient if the rotation profile and the 
shear-field strength in the source region can be estimated (see Appendix B). 

viii) The relative contributions to the total rotation by the shear field and by the twist can 

be disentangled by comparing both the observed rotation and rise profiles with the 
corresponding curves from a model, since these profiles possess a different dependence 
upon the parameter combination. The resulting estimate for the twist allows one 

to judge the occurrence of the helical kink instability. 

ix) Magnetic reconnection contributes only weakly (much less than the shear field) to the 
total rotation in the simple bipolar source regions considered. 

From the comparison with the simulation of rotating flux in Lynch et al. (2009) we con- 
clude: 

x) The rotation profile differs strongly between configurations with and without a flux rope 
at the onset of the eruption. 

The comparison with the stereoscopic observations and three-dimensional reconstruction 
of the erupting prominence in the 9 April 2008 “Cartwheel CME” additionally shows the 
following. 

xi) The rotation profile obtained in Paper I from the stereoscopic reconstruction of STEREO 
data is equally well reproduced by our model up to heights & 1.5 R Q above the photo- 
sphere for a range of d>-B e t combinations which include a strongly kink-unstable case 
(cf> = 57 t, B et /B ep = 0.42), a weakly kink-unstable case (<t> = 3.5;r, B et /B ep = 0.67), 
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and a kink-stable case (<£> = 2.5tt, B et /B ep = 1.06). However, the strongly kink-unstable 
configuration is ruled out by the simultaneous consideration of the rise profile, and sev- 
eral features of the kink-stable model argue strongly against this configuration. These 
are the implied high value of the shear field, the rotation profile at greater heights, and 
the unrealistic start height of the unstable rise of ~ 0.2 R Q . Hence, the occurrence of a 
weak helical kink instability in the Cartwheel event is very likely. 

Our results add to the complexity of the phenomenon of flux rope rotation in eruptions 
which is already known from investigations that focused on the influence of reconnection 
(e.g., Jacobs et al., 2009; Shiota et al, 2010; Cohen et al, 2010; Thompson, 2011; Lugaz 
et al., 2011). An overall very complicated dependence on several parameters and on the 
structure of the ambient field is revealed. Thus, the quantitative prediction of the rotation 
is a difficult task. The parametric study performed here indicates for simple bipolar source 
regions that the strength of the external shear field is the primary parameter determining the 
total rotation. The flux rope twist and the height profile of the external poloidal field are 
of relatively minor importance as long as they stay in the typical ranges indicated by the 
observations. We did not yet study a possible influence of the height profile of the external 
shear field. The external shear field of filament channels may be estimated to sufficient 
precision from a simple linear force-free field extrapolation. It will be worth testing whether 
numerical modeling starting from such fields, embedded in current-free outer field, yields 
rotations in agreement with observations of eruptions from bipolar source regions. 

Several investigations indicate that erupting flux ropes align with the heliospheric current 
sheet in the course of their interplanetary propagation (e.g., Bothmer and Schwenn, 1998; 
Yurchyshyn, 2008; Paper I). This suggests that the coronal rotation merely decides whether 
a parallel or an antiparallel alignment will result at 1 AU. However, since complex physics 
is involved and since rotations on the order of 90° may not be rare, the quantitative study of 
the effects that determine the rotation in the corona remains of high scientific and practical 
interest. 

Acknowledgements We acknowledge the careful reading of the manuscript and the constructive com- 
ments by the referee. BK acknowledges support by the DFG, the STFC, and by NASA through Grant 
NNX08AG44G. TT’s work was partially supported by the European Commission through the SOTERIA 
Network (EU FP7 Space Science Project No. 218816) and by the NASA HTP and LWS programs. WTT’s 
work was supported by NASA Grant NNG06EB68C. 


Appendix A: Influence of the Coronal Polarity Inversion Line 

There are quite strong indications that CMEs align with the heliospheric current sheet in the 
course of their propagation, i.e., with the PIL in the solar wind (see references in Section 5). 
This leads to the question whether the PIL guides the rotation of erupting flux ropes also in 
the corona. Here the PIL formed by the external field, due to sources outside the flux rope, 
must be considered. We use “CPIL” to denote this structure in the corona, where [3 < 1. The 
heliospheric current sheet and the CPIL differ in two properties of relevance here. Lirst, in 
the solar wind p > 1 , so that the pressure gradient is generally dominant over the Lorentz 
force, while the opposite is true in the corona. Second, the heliospheric current sheet is 
the location of pressure gradients and Lorentz forces, while the CPIL generally lacks both. 
In the low-beta corona, currents are induced at separatrix surfaces, or at quasi- separatrix 
layers, if the equilibrium is perturbed or lost. The CPIL generally does not coincide with 
these structures. Therefore, the CPIL should not influence the rotation of erupting flux ropes 
in this height range. 
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Figure 10 Visualization of the weakly kink-unstable modified Titov-Demoulin equilibrium (<J> = 3.577, 
Bet / B e p = 0.67; Figures 2 and 4) whose eruption characteristics match the observations of the Cartwheel 
CME best (top panel) and of the corresponding external field (bottom left) and potential field (bottom right). 
The magnetogram and field lines starting in the photospheric flux concentrations are shown. 


Figure 11 Orientation of the 
PIL in the external field of the 
configuration shown in Figure 1 0 
at the position of the flux rope 
and different heights. The 
orientation is indicated by a black 
line. 



Figure 10 shows the initial equilibrium of the weakly kink- unstable run which matches 
the Cartwheel event best, the corresponding external field, and the potential field that results 
when the full magnetogram of the vertical field component of the equilibrium, B : (x, y, 0), is 
extrapolated into the volume above. The full magnetogram includes the contributions from 
the flux rope, which are excluded from the external field. The CPIL of this configuration 
at the photospheric and three coronal levels is shown in Figure 1 1 . The CPIL changes its 
orientation in a clockwise sense if one goes upward, but the unstable flux rope rotates in a 
counterclockwise direction, since it is left handed. The clockwise changing CPIL orientation 
results from the dominance of the external toroidal field, B et , over the external poloidal 
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Figure 12 Same as Figure 1 1 
for the potential field of the 
configuration shown in Figure 10. 



field, B ep , at great heights. This situation can typically be expected to occur because B cl 
typically has a larger spatial scale than B ep (set by the distance between the sources in 
the photosphere). The important fact here is that the CPIL does not appear to have any 
significant influence on the rotation of the flux rope in the zero-beta simulations performed 
in this paper. For the reasons given above, this is valid also if other height profiles of B cx or 
B e p lead to a different profile of the CPIL orientation with height. 

We also consider the approximation of the true CPIL by the PIL in a potential-field 
extrapolation of the full photospheric magnetogram, B z (x, y, 0). In practice, it is difficult 
or even impossible to determine the external field. This requires the determination of the 
coronal currents through a nonlinear force-free extrapolation from a vector magnetogram. 
The former is still difficult to carry out and the latter may not be available. The PIL in the 
potential field extrapolated from the magnetogram of the weakly kink-unstable configuration 
in Figure 10 is shown in Figure 12. Its orientation vs. height is very similar to the behavior of 
the true CPIL. This supports the conclusions drawn in Paper I from a potential- field source- 
surface extrapolation for the source region of the Cartwheel CME. 


Appendix B: Inferring the Helical Kink Instability 

Finally, we consider the question how the occurrence of the helical kink instability can be 
inferred in spite of the typical dominance of shear-field-driven rotation above Kl-driven 
rotation in CMEs, which is suggested by the present study. The writhing of erupting flux 
loops into a helical shape, i.e., the apex rotation, has often been regarded to be a clear 
indication if not a proof of the occurrence, but from Isenberg and Forbes (2007), Lynch 
et al. (2009), and the present investigation it is obvious that the apex rotation, by itself, 
is not conclusive in general. To infer the occurrence unambiguously, one needs to find a 
supercritical amount of twist or features incompatible with the shear-field-driven writhing 
(if other rotation processes can be excluded, such as reconnection with the ambient field and 
the straightening from a strong initial S shape). 
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Figure 13 Flux rope rotation vs. 
normalized apex height for a 
range of Bet-® combinations 
that yield stable flux rope 
equilibria near the threshold of 
the helical kink mode. The 
normalized sunspot 
semi-distance is L = 0.83 as in 
Figure 6. Observed rotation 
profiles which exceed the 
rotations plotted here for the 
same ratio B e t/B e p indicate the 
occurrence of the helical kink. 



The twist can be calculated from extrapolated magnetic fields, or be estimated from the 
observation of twisted substructures or from parametric studies like the one performed here. 
The second option yields indications for supercritical values only rarely {e.g., Vrsnak, Ruzd- 
jak, and Rompolt, 1991; Romano, Contarino, and Zuccarello, 2003), and the other options 
are computationally demanding. In particular, the extrapolation must at least be based on 
the nonlinear force-free field assumption (Leka, Fan, and Bames, 2005; Valori et al., 2010), 
which additionally requires a vector magnetogram. 

Features incompatible with purely shear-field-driven rotation are approaching flux rope 
legs, an apex rotation considerably exceeding 90°, and the development of more than one 
helical turn in the shape of the erupting structure as a whole (i.e., in its axis). The shear 
field causes the flux rope legs to lean aside in opposite direction perpendicular to its own 
direction, i.e., to move away from each other. The legs of a kinking flux rope approach each 
other if the initial twist exceeds ~ 6n (Kliem et al, 2010). However, this will occur only 
rarely, due to the high twist required. 

The rotation resulting from the opposite bending of the flux rope legs by the shear field 
would always stay below 90° if there were no further contributing effect. As discussed in 
Section 3.1, an additional contribution from the relaxation of the twist in a force-free flux 
rope does, in fact, always exist, even if the twist is insufficient to trigger the helical kink. This 
contribution reaches ~ 40° if the shear field is very small (Figure 6) (and the sunspot distance 
stays within the typical range observed; see Section 3.2 and Figure 7). It is expected to 
remain smaller when a larger shear field reduces the twist through its own enforced rotation. 
Consequently, a rotation exceeding 1 30° (in the absence of other contributions) is a strong 
indication that the helical kink instability occurred. 

Figure 6 also shows that the instability is indicated already by somewhat smaller values 
of the rotation if the angle rises quickly at low heights. To quantify this, we have computed 
the family of rotation curves shown in Figure 13, which all represent initial conditions stable 
against the helical kink mode but close to its marginal stability boundary in the B e t -d> plane. 
Any rotation profile lying clearly above the corresponding curve in this figure, with the shear 
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field set to the value in the eruption’s source region, indicates a contribution by the helical 
kink. This requires estimating the shear field B et and the footpoint distance, 2D{, of the 
erupting flux. A linear force-free field extrapolation yields a reliable estimate of the shear 
field, and the consideration of the distances of the main polarities along and across the PIL 
may yield a rough estimate in source regions of relatively simple structure. 

These considerations are valid only as long as no other processes contribute to the rota- 
tion, which could be reconnection, the straightening from a strong initial S shape, and, at 
heights above the corona, the alignment with the heliospheric current sheet. The straighten- 
ing from an initial S shape can simply be subtracted from the total rotation angle before the 
comparison with the plots in Figure 13 is made. 

Since the shear-field-driven writhing is caused by antiparallel forces on the legs, the 
erupting flux rope’s axis will never develop more than one helical turn. Some eruptions 
indicate more than one turn, from an S shape of one or both legs {e.g., Ji et al., 2003; 
Maricic et al., 2004) or from a very writhed shape as a whole {e.g., Romano, Contarino, and 
Zuccarello, 2003). The helical kink instability can be inferred in such cases from the shape 
of the structure (Torok and Kliem, 2005; Kumar et al, 2012). 
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